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I. INTRODUCTION 



The 2D g-state Potts models for various q have been of interest as examples of 
different universality classes for phase transitions and, for q = 3, 4, as models for the adsorp- 
tion of gases on certain substrates. The q = 2 Ising special case has long served as a simple 
exactly solvable model of cooperative phenomena. However, for q > 3, the free energy of 
the Potts model has never been calculated in closed form for arbitrary temperature. It is 
thus worthwhile to obtain further information about the properties of the Potts model, and 
we shall do this in the present paper via calculations of complex-temperature (CT) zeros of 
the partition function of the Potts model for the honeycomb and kagome lattices. One of 
the motivations for this work is the recent calculation and analysis of long low-temperature 
series for the g-state Potts model on these lattices by Jensen, Guttmann, and Enting 
Our results enable one to relate the CT singularities in thermodynamic quantities found in 
Ref. |§ to positions on the CT phase boundaries of the respective models. 

The study of statistical mechanical models with magnetic field ||, temperature ||- [|J, 
or both H generalized from real to complex values has yielded interesting insights into 
the properties of these models. For a discrete spin model at temperature T and in an 
external magnetic field H, the partition function Z is, up to a prefactor, a polynomial in 
the Boltzmann weights z(K) and fi(h) containing dependence on K = (3J and h = /3H, 
where (3 = 1/(/cbT), and J is the spin-spin coupling. It is of interest to study the zeros 
of Z (i) in the complex \i plane for physical T H; (ii) in the complex z plane for physical 
(vanishing or nonvanishing) H ||; and (iii) on the C 2 manifold (/x, z) when both K and h are 
complex ||. Here we shall concentrate on case (ii), i.e., Fisher zeros. In the thermodynamic 
limit, via a coalescence of zeros, there forms a continuous locus B of points where the free 
energy is nonanalytic. This locus serves as the union of boundaries (whence the symbol 
B) of the various complex-temperature phases. Thus, calculations of complex-temperature 
partition function zeros on sufficiently large finite lattices yield useful information on the 
CT phase diagram in the thermodynamic limit. (Hereafter, to avoid repetition, we shall 
simply refer to zeros of the partition function, it being understood that these are complex- 
temperature zeros.) In making inferences from such finite- lattice calculations about B in the 
thermodynamic limit it is important to vary both the lattice size and the type of boundary 
conditions to have an accurate idea of the sensitivity of the locations of the zeros to these 
choices. Some of the earliest work on CT properties of spin models dealt with these zeros [^,§J. 
Another major reason for early interest in these properties of spin models was the fact that 
unphysical, CT singularities complicated the analysis of low-temperature series expansions to 
get information about the location and exponents of the physical phase transition || . A third 
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reason for interest in these properties is the fact that, as additional sources of information 
about thermodynamic functions, they can expedite progress toward exact solutions. Aside 
from well-understood exceptions [pLOfl , CT singularities of thermodynamic functions occur 
on the continuous locus of points B where the free energy is nonanalytic. Hence, when 
investigating these singularities, it is useful to do so in conjunction with a calculation of the 
zeros of the partition function to infer the approximate location of the phase boundary B 
separating various CT phases [Q. Interestingly, some of these singularities can be related 
directly to physical singularities: by using duality, one can show an exact equivalence of 
the free energy of the g-state Potts antiferromagnet on a lattice A for the full temperature 
interval < T < oo and the free energy of the g-state Potts model on the dual lattice for a 



semi-infinite interval of complex temperatures ||13|| . This implies the existence of two quite 
different types of CT singularities: the generic kind, which does not obey universality or 
various scaling relations ||T4| , p!5| , and a special kind which does obey such properties and 
encodes information of direct physical relevance. Although we consider the honeycomb and 
kagome lattices here, we mention that previous calculations of zeros of the partition function 



for the Potts model with q > 3 have been done on the triangular and square lattices |[16Hl9 



II. MODEL 

The (isotropic, nearest-neighbor) g-state Potts model on a lattice A is defined by the 
partition function 

with the Hamiltonian 

(nn 1 ) n 

where a n = l,...,q are Z g -valued variables on each site n G A, (3 = (/c^T) -1 , and {nn') 
denotes pairs of nearest-neighbor sites. The symmetry group of the Potts Hamiltonian is the 
symmetric group on q objects, S q . We use the notation introduced above, K = (5 J, h = (3H, 
and 



a = z- 1 = e K 



(2.3) 



The (reduced) free energy per site is denoted as / = —f3F = lini7v s -^oo A^T 1 In Z, where N s 
denotes the number of sites in the lattice. There are actually q types of external fields which 
one may define, favoring the respective values a n = 1, ..,(?; it suffices for our purposes to 
include only one. The order parameter (magnetization) is defined to be m = (qM—l)/(q—l). 
where M = (a) = lim^o df/dh. With this definition, m = in the symmetric, disordered 
phase, and m = 1 in the limit of saturated ferromagnetic (FM) long-range order. We 
consider the zero-field model, H = 0. For J > and the dimensionality of interest here, 
d = 2, the g-state Potts model has a phase transition from the symmetric, high-temperature 
paramagnetic (PM) phase to a low-temperature phase involving spontaneous breaking of the 
S q symmetry and onset of ferromagnetic (FM) long-range order. This transition is continuous 
for 2 < q < 4 and first order for q > 5. The critical exponents and universality classes of the 
cases where the model has second-order transitions are well understood P,|20|. The g-state 
Potts model has the property of duality |I|,|2] j2l|j22]| , which relates the partition function on 



a lattice A with temperature parameter a to the partition function on the dual lattice with 
temperature parameter 

a + g — 1 1 

a d = V(a) = — , i.e. x d = - (2.5) 

a — 1 x 

Other exact results include formulas for the PM-FM transition temperature on the square, 
triangular, and honeycomb lattices |I|j2]j22]| , and calculations of the free energy at the phase 



transition temperature, and of the related latent heat for q > 5 [23|. No formula is known 



for the PM-FM transition on the kagome lattice, although there have been a number of con- 
jectures; for a recent discussion, see Ref. 0. The case J < 0, i.e., the Potts antiferromagnet 
(AF) has also been of interest because of its connection with graph colorings and the fact 
that, for certain lattices and values of q, it exhibits nonzero ground state entropy [24 , 25] ; for 
a recent discussion, see [0,0 an d references therein. Depending on the type of lattice and 
the value of q, the model can also have a phase with AFM long-range order. For q > 3 on the 



honeycomb lattice there is no AFM phase p8| , p9| . For any lattice A, the partition function 
can be expressed in a form involving a sum of powers of q which allows a generalization from 
positive integer q to real (or, indeed, complex) q, and we shall use the generalization to real 
q at certain places below. Reviews of the model include Refs. HJTTjI. 

On a finite lattice, the g-state Potts model partition function Z is a polynomial in the 
Boltzmann weight a. We calculate this polynomial by transfer matrix methods. This is a 
challenging numerical problem for large lattices, since the degree of the polynomial is equal 
to the number of bonds, = (A/2)N S , where A is the coordination number, and there is 
a very large range in the sizes of the coefficients, from q for the highest-degree term a Nb to 
exponentially large values for intermediate terms. The latter property is obvious from the 
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fact that for K = 0, i.e., a = 1, the sum of the coefficients in Z is q Ns . From this, we then 
compute the zeros. A general property of the CT phase boundary for any lattice and q value 
is invariance under complex-conjugation: B — ■> B as a — > a*. 

In addition to the locations of the curves comprising the CT phase boundary B inferred in 
the thermodynamic limit from the zeros calculated on finite lattices, one can extract further 
information. As one approaches the thermodynamic limit, so that one can define a density 
of zeros, this density normally behaves near a singular point a s as 

g{s) ~ s 1- " , as s -> (2.6) 

where s denotes the arclength along B away from a s (so that s — \a — a s \ as s — > 0) and 
where the singularity in the free energy at a s is / S j n5 ~ \a — a s | 2 ~ a If the partition 

function has a zero at some point zq with a multiplicity proportional to the number of lattice 
sites, N s , then this formula, eq. (2Ji), is modified by the addition of a term proportional to 



a delta function 5(s). In Ref. [FI it was proved (as Theorem 6) that for the Ising model on 
a lattice with odd coordination number, this happens at z — — 1. In particular, this occurs 
for the Ising model on the honeycomb lattice (see further below). 



III. PARTITION FUNCTION ZEROS ON THE HONEYCOMB LATTICE 

A. Comparison with Exact B for Ising Case 

In order to study the effects of the finite lattice size and of different boundary conditions, 
as well as checking the computer programs used, it is valuable to calculate the zeros for the 
q = 2 Ising case where the resulting locus of zeros can be compared with the exactly known 
CT phase boundary B. As noted above, these zeros, like the others to be presented further 
below, are calculated by a transfer matrix method. From the known expression for the free 
energy, this boundary was determined in Ref. it is the locus of solutions to the equation 

1 - 2a + 6a 2 - 2a 3 + a 4 - 2a(l - a) 2 p = (3.1) 

where —3/2 < p < 3 f32fl . This locus is shown in Figs. H and Because q = 2 and 
the honeycomb lattice is bipartite, the CT phase boundary B and also the set of zeros are 



invariant under the inversion map [33]. The CT boundary consists of the union of two parts. 
The first is an arc of the unit circle extending from 9 = arg(a) = tt/3 around through a = — 1 
to 9 = — 7r/3, while the second is a lima bean-shaped curve that crosses the positive real axis 
at the PM-FM critical point, a PM -FM, q =2 = 2 + ^3 = 3.732... and at the PM-AFM critical 
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point, apM-AFM, q =2 = a pM-FM,q=2 = % ~ v^3 = 0.267949.... These two parts intersect each 
other at multiple points at ±z; these multiple points are singular points of B in the sense 
of algebraic geometry [35]. The phase surrounding the origin in the a plane is the AFM 
phase; the one surrounding the infinite-temperature point a = 1 is the PM phase, and the 
one extending outside B to complex infinity in all directions is the FM phase [12| . 



Before we start to present our results, we have to introduce our notation for the sizes 
and orientations of the lattices. We recall first that the (infinite) honeycomb lattice is a 
homopolygonal member, and the kagome lattice, a heteropolygonal, member of the class 
of Archimedean lattices, i.e., regular tiling of the plane by one or more types of regular 



polygons such that every vertex is equivalent to every other vertex [36|. An Archimedean 



lattice is thus uniquely defined by the ordered sequence of polygons that one traverses in 
making a circuit of any vertex. In standard mathematical notation [36|, such a lattice is 
denoted A = (n^D where Pi refers to the type of polygon and a« denotes the number 
of times that it appears consecutively in the product. In this notation, the kagome lattice 
is denoted (3 ■ 6 • 3 • 6). The homopolygonal subset of Archimedean lattices (consisting of 
tilings with only one type of regular polygon) is closed under duality, but the heteropolygonal 
Archimedean lattices (consisting of regular tilings using more than one type of polygon) have 
duals that are not Archimedean lattices. In particular, the dual of the kagome lattice, called 
the diced lattice, is not Archimedean. This mathematical background will be useful below 
when we give results for the diced lattice. To indicate the size of a given lattice for both the 
honeycomb and kagome cases, we count the number of hexagons. As an illustration, the sizes 
of the honeycomb and kagome lattices in Fig. [I] are 4x3 and 3x4 hexagons, respectively. 
The number of sites in a lattice is also dependent on the boundary conditions: with periodic 
boundary conditions in the horizontal direction for example, the sites on the left and right 
are identified, while with free boundary conditions they are counted independently from each 
other. As is evident from Fig. [I], a honeycomb lattice of size, in our notation, N x x N y , is 
maximally square-like if one takes N x slightly larger than N y . 



FIG. 1. Honeycomb and kagome lattices to illustrate our conventions for indicating sizes. 
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Since we use duality at many points in this work, we chose lattices that have natural 
dual lattices. This excludes lattices that are periodic in both directions, for the following 
reason: duality relies on the fact that every closed polygon divides the lattice into at least 
two regions. However, a lattice with periodic boundary conditions in both directions, and 
hence with toroidal geometry, has the property that there exist closed contours that do 
not divide the surface into two disjunct regions. Since boundary effects are, in general, best 
suppressed if one uses periodic boundary conditions in as many directions as possible, we use 
boundary conditions that are periodic in one direction and free in the other, i.e., cylindrical 
boundary conditions. Our notation for the boundary conditions (BC's) is (fbc,pbc) for free 
and periodic BC's in the horizontal (x) and vertical (y) directions, respectively (see Fig. |l|), 
and (pbc,fbc) for periodic and free BC's in the x and y directions. This notation makes 
explicit the direction in which the cylindrical boundary conditions are periodic. 

Complex-temperature partition function zeros for the q = 2 Ising case of the Potts model 
on the honeycomb lattice are shown in Fig. ^| for both (fbc,pbc) and (pbc, fbc). The gray 
curves are the exactly known CT phase boundary B. We remark on several features. 

• The partition function has a multiple zero at z — — 1 with multiplicity oc N s for large 
N s This follows from Theorem 6 of Ref. |TTJ and corresponds to the term oc (1 — z 2 )~ 2 



in the expression for the specific heat C in the FM phase given as eq. (3.12) in Ref. 
PH ; the apparent additional singularity at the infinite-temperature point z — 1 is not 
relevant since the formula does not apply in that region. 

The zeros lie very close to the arcs protruding into the PM phase. 

There seems to be some repulsion of the zeros from the multiple points at a = ±i 
(similar to what was seen in Ref. ||). 

For the (fbc, pbc) case, the zeros on the unit circle show no radial deviation. 

In general, the zeros calculated with the choice (fbc, pbc) lie closer to the exact bound- 
ary curves B than those calculated with (pbc, fbc). 

The zeros lie on the outer side of the boundary between the PM and FM phase. 
This can be understood as a consequence of the fact that with either the (fbc, pbc) 
or (pbc, fbc) boundary conditions, the sites on the free boundary have a coordination 
number of two rather than the usual A = 3 for sites on an infinite honeycomb lattice. 
Hence, ordering effect of the spin-spin interactions is commensurately reduced, thereby 
reducing the finite-lattice manifestations of the ordered, FM phase, i.e. shifting the 
PM-FM boundary outwards. 
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FIG. 2. Partition function zeros in the complex a plane for the Ising (q = 2 Potts) model on a 

honeycomb lattice, (a) left: size 9 x 12 hexagons and (fbc, pbc); (b) right: 10 x 10 hexagons and 

(pbc, fbc). 

Note that for the (fbc, pbc) choice, there is one site per hexagon at the boundary with only 
two instead of the usual A = 3 bonds. For the (pbc, fbc) BC's, there are two of these sites 
for each of the hexagons on the upper and lower boundaries. This motivated us to formulate 
a third kind of boundary condition: starting from the (pbc, fbc) choice, we added bonds 
connecting the boundary sites with fewer than three bounds so that all sites on the lattice 
have the same coordination number A = 3. We denote this choice as (pbc,fbc)A- The zeros 
calculated with this third choice of boundary conditions are plotted in Fig. |3|. The main 
difference relative to the previous two choices of BC's is that the zeros in the Re(a) > 
half plane have less scatter, lie closer to the exact boundary B, and also, in some cases, lie 
inside the PM-FM phase boundary. The zeros with Re(a) < are less scattered than those 
with the (pbc, fbc) choice and track the arc of the unit circle well, although they do not, 
in general, lie on it, as was the case with the choice (fbc,pbc). The conclusion from this 
comparison with exactly known results is that, if one did not know the exact boundary B 
to begin with, one would be able, by combining results on zeros calculated with different 
boundary conditions, to reconstruct it with reasonable accuracy. 

The density of zeros on B near the physical PM-FM transition is consistent with vanishing 
according to eq. (|2.6f) with a = 0, i.e., g ~ \a— apM-FM,q=2\ as \a — apM-FM,q=2\ — > 0. By the 
a — > 1/a symmetry, the same is true of the PM-AFM transition. The situation at z = —1 is 
more complicated because the partition function has an isolated zero of multiplicity scaling 
like the lattice size there; consequently, as discussed above, the density g has a delta function 



7 



term oc 5(s) as well as its usual term (|2.6|) , where s denotes the arclength on B away from 
the point z = —1. The analysis of Ref. |31[] found that at z = —1, aside from the leading 
singularity ~ (l + z)~ 2 in the specific heat, there is also a subleading logarithmic divergence; 
it follows that the density of zeros on B near to z = — 1 has, in addition to the delta function 
term, a term that vanishes like s. The zeros in Figs. §-|3| are consistent with this. 
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FIG. 3. Partition function zeros for the Ising model on the honeycomb lattice of size 10 x 10 
hexagons and (pbc, fbc)A BC's. 




B. q = 3 Case 

For general q, from duality and a star-triangle relation, an equation yielding the value of 



the PM-FM transition point has been derived p2j , viz., x —3x — y/q = 0, or, in terms of a 



a 3 -3a 2 -3(g-l)a-g 2 + 3g-l = (3.2) 

For < q < 4, this equation has three real roots, while for q > 4 (and the formal values 
q < 0) it has one real root. The motion of the real roots as a function of q is plotted in Fig. 
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FIG. 4. Real roots of eq. Q3,2| ), as a function of q. 
For q = 3, the solutions are 

ai,a=3 = flpM-FM,= 3 = 1 + 2^008(^/18) = 4.41147... (3.3) 



02,9=3 = 1 - v 7 ^ cos(tt/18) + 3 sin(vr/18) = -0.1847925... (3.4) 



and 



a 3q=3 = 1 - v 7 ^ cos(tt/18) - 3 sin(7r/18) = -1.22668... (3.5) 
The point a PM _ FM;q=3 is the physical PM-FM critical point. As discussed in Ref. [E9J, if one 



follows the roots of eq. ( |3.2| ) as q is changed continuously, one sees that the middle root a 2 
decreases from the PM-AFM critical point 2 — \f?> for q = 2 through at q z = (3 + v^5)/2 = 
2.618... to the negative value Q3.4p for q = 3. This reflects the fact that as q increases from 
2 to q z , the physical AFM phase is squeezed out. 

Our zeros of the partition function for the q = 3 case on the honeycomb lattice are shown 
in Figs. H - ^ for the three types of boundary conditions discussed before. 
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FIG. 5. Partition function zeros in the a plane for the q = 3 Potts model on a honeycomb lattice 
of size 8x6 hexagons and boundary conditions of type (fbc,pbc). 
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FIG. 6. Partition function zeros in the a plane for the q = 3 Potts model on a honeycomb lattice 



of size 8x6 hexagons and boundary conditions of type (pbc,fbc). 
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FIG. 7. Partition function zeros in the a plane for the q = 3 Potts model on a honeycomb lattice 
of size 8x6 hexagons and boundary conditions of type (pbc,fbc)A- 

In all three plots, the roots apM-FM, q =3 and a 3tq=3 correspond very well to points where 
the CT phase boundaries comprising B (as inferred for the thermodynamic limit from these 
zeros on finite lattices) cross the real axis. Hence, one anticipates that a CT phase boundary 
might cross the real axis at the value of the middle root, 02,5=3 = —0.1848.... From the zeros 
in Figs. |5| - [7], one can indeed infer that in the thermodynamic limit a CT phase boundary 
curve could cross the real axis at this point. Since the specific heat exponent has the known 
value a = 1/3 for this model 0, it follows from eq. ( |2.6| ) that the density of zeros near the 
physical PM-FM transition point vanishes like g(s) ~ s 2//3 , where s is the arclength on B 
away from this point. This is consistent with the calculated zeros. 

Further, we see arcs protruding into the PM phase, ending at complex conjugate (c.c.) 
points a e ,a* = 0.37(2) ± 1.29(3)i, where the numbers in parentheses refer to the estimated 
uncertainties in the final digits. Evidently, these are the analogues for q = 3 of the exactly 
known arcs for the Ising q = 2 case. While the arcs in the Ising case have endpoints on the 
unit circle at angles 9 = ±7r/3, the endpoints in the q = 3 case lie slightly farther out from 
the origin, at \a e \ ~ 1.3, and have slightly larger angles 9 ~ ±75°. We find that this trend 
is true for larger q values also, i.e., \a e \ and arg(a e ) increase with increasing q. 

In addition, there are at least two more points at which curves of zeros cross the real a 
axis, at a = = —2.77(3) and at a = —0.65(2). In Ref. |13| it was shown that if the g-state 
Potts antiferromagnet on the dual lattice has a PM-AFM transition at a p M—AFM,k d 1 then 
the dual image of this, namely, T>(a PM _AFM,A d ) = cli,a, is the leftmost point in the a plane 
where B crosses the real axis. Since the PM-AFM point satisfies < a p m—af M,k d < 1, it 
follows that the dual image gi^a satisfies —00 < ci^a < —(q— 1). In particular, for q = 3, this 
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connection was used, in conjunction with a precisely measured value of a, p m—af M,t,q=3 on the 



triangular (t) lattice [37] to infer the value of ag for the model on the honeycomb (he) lattice: 
ahc,e, q =3 = T>(ap M -AFM,t,q=3) = ~ (2.76454 ± 0.00015) , i.e., z hcAq=3 = a^ Aq=3 = -(0.36172 ± 
0.00002). Our zeros are in agreement with this result. This point also manifests itself as a 
singularity evident from low-temperature series for the specific heat C, magnetization m, and 
susceptibility Xi which yield the value z = —0.363 ± 0.003 f3j. Using duality and the weakly 
first order nature of the physical PM-AFM transition of the q = 3 Potts antiferromagnet, 
it follows that the free energy also has the same weakly first order singularity at a hc ^ q=3 . 
The low-temperature series analysis of Ref. found evidence for a continuous transition at 
this point, with exponents ag = 0.5, = 0.11, and 7^ = 1.15. We have repeated the series 
analysis with dlog Pade approximants (PA's) and differential approximants (DA's) []38[| . Our 



DA results also yield a ~ 0.5; our PA's did not locate the singularity with sufficient precision 
to infer a reliable value for a. Given the duality and the fact that a = a 1 for the physical 
PM-AFM transition of the q = 3 Potts AF on the triangular lattice, it follows that the 
singularity in the free energy of the q = 3 Potts model on the honeycomb lattice at must 
also be the same as approached from the right or left. Since the singularity in the internal 
energy at a singular point a s is U S i ng ~ \a — a s | 1_a , one normally assigns the formal value 
a = 1 to a first-order transition. A possible way of reconciling these results is to observe 
that if a first order transition occurs superimposed with a divergent specific heat, then one 
could get a value of a < 1 in fitting the transition. For example, consider an illustrative 
internal energy function that behaves near a phase transition point like 

U(T) ~ U analytlc + c 1>+ 0(T - T c ) + c 2 , + |T - T c \ 1 ' 2 (3.6) 

for T \ T c , and similarly for T /* T c , with the coefficients replaced by Ci ( _ and C2,_. Here, 
U analytic denotes terms that are analytic near T c and Q(x) is the step function, 6(x) = 1 
if x > and otherwise. As one approached T c from above (below) a high-temperature 
(low-temperature) series analysis would give a = 1/2, but the transition would still be first 
order because of the discontinuous term. A one-sided version of this behavior occurs in the in 
the six-vertex model for the ferroelectric compound potassium dihydrogen phosphate (KDP) 
|39| ; in that case, the form (|3.6|) applies for the high-temperature side, while U is a constant 
on the low-temperature side. 

Another source of information on a is the density of zeros. However, it is difficult to 



use this to obtain an accurate value of a. For example, Ref. [1£] included calculations of 
zeros for the g-state square-lattice Potts model not just for the values q = 3,4 where the 
PM-FM transition is continuous, but also for the values q = 5 and 6, where this PM-FM 
transition is first order; see Figs. 3 and 4 therein. For these cases, one would formally set 
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a = 1 as mentioned above, so that eq. (|2.6| ) would predict that the density g of zeros should 
remain essentially constant up to the endpoint of the distribution (of course, the positivity 
of the coefficients of the terms in the partition function means that for a finite lattice, there 
cannot be any zeros on the positive real axis in the a or z plane). This is consistent with 
the plots of zeros for these q = 5 and 6 cases, but it would be difficult to extract an accurate 
estimate of a from those plots. Below we shall present a similar plot for another case where 
the PM-FM transition is known to be first order, namely the q = 5 Potts model on the 
honeycomb lattice, and a similar comment applies to this plot. 

Further CT singularities and their relation with the boundary B will be discussed else- 
where in work with the authors of Ref. 0. Our calculations also suggest that there are 
several unphysical O phases which overlap with parts of the negative real a axis. There may 
be other O phases that do not touch the real axis, but the resolution is not high enough to 
make a definitive statement here. 

Concerning the sensitivity of the zeros to lattice boundary conditions, several remarks are 
in order. The zeros in the Re(a) > half plane are relatively insensitive to these boundary 
conditions. However, certain features of the zeros in the Re(a) < half plane do show such 
sensitivity. This is similar to what was found from a comparative study of different boundary 
conditions for the zeros of Z for the g-state Potts model on the square lattice for several 
values of q |7|,[19| (see also Ref. @). 

We have carried out similar calculations of zeros for the q = 4 Potts model on the 
honeycomb lattice, and these will be reported in joint work with the authors of Ref. ||. 
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C. q = 5 Case 



Ira(a) 0- 



Re(a) 

FIG. 8. Zeros of Z in the a plane for the q = 5 Potts model on a honeycomb lattice of size 7x6 
hexagons and boundary conditions of type (fbc,pbc). 



It is also of interest to investigate a value of q in the range where the PM-FM transition 
is first order, i.e., q > 5. We have done this for the value q = 5, and we show a resulting 
plot of zeros in Fig. || Here, eq. ( |3.2| ) has the single real root, which is the PM-FM critical 
point, 



®PM-FM,q=5 



2-1/3 . 5 1/2 (1 + 51/2)1/3 + 2 l/3 . 5 1/2 (1 + 5 l/ 2) -l/3 + l = 5.5398... (3.7) 



Since the Potts antiferromagnet with q = 3 and q = 4 on the triangular lattice has, respec- 
tively, a finite-temperature PM-AFM phase transition [0,0 anc ^ a zero-temperature critical 
point P5[ , it is expected that for q > 5, the model is disordered for all temperatures. This, 
together with the connection discussed in Ref. [jTj|, would imply that the leftmost point at 
which B crosses the real axis for the q = 5 Potts model on the honeycomb lattice is ag < —4. 
Our zeros are consistent with this. 



D. Further Discussion 

It is a general feature of the maximal (or sole) real solution of eq. ( ^2|) , i.e., cipm-fm, 
that it increases monotonically with q for q > 0. This is evident in Fig. § and reflects 
the basic thermodynamic property that as q increases, the spins become "floppier" , and one 
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must go to lower temperature to obtain FM long range order. In addition to the features 
already discussed, we note that (i) the leftmost point where B crosses the real a axis, a,£, 
moves to the left as q increases; and (ii) the points where B crosses the imaginary axis move 
out from the origin as q increases. Both of these features can be understood, as discussed 
directly above, by the reduction in the size of the (FM) ordered phase as q increases. 

Because of the duality relation, these partition function zeros, in the a plane, of the re- 
state Potts model on the honeycomb lattice also yield equivalent zeros of the same model on 
the dual, triangular lattice in the plane of the variable given in eq. ( |2.5|) . A comparison of 
the plots calculated with different boundary conditions is valuable since this gives a measure 
of the effects of these boundary conditions (see also Refs. |T6|j40[|). 



IV. PARTITION FUNCTION ZEROS ON THE KAGOME LATTICE 



A. Comparison with Exact B for Ising q = 2 Case 
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-2-10 1 2 

Re (a) 

FIG. 10. Zeros of Z in the a plane for the Ising model on a kagome lattice with 4x6 hexagons 
and (pbc,fbc)A boundary conditions. 

For our calculations of zeros of Potts model on the kagome lattice, two of the boundary 
conditions that we use are (pbc,fbc) and (pbc,fbc)A, where now A = 4. The third can be 
described as follows: we start on the lattice that is the dual to kagome, namely the diced 
lattice, and impose (pbc,fbc) boundary conditions; then we transform the results back to 
the kagome lattice by the duality map on a, eq. ( |2.5| ). To save space, for each value of q, 
we only show results for the first two of these choices of boundary conditions. For q = 2 



these are given in Figs. |] and [lCj. The exact CT phase boundary B is given by the locus of 
solutions of the equation 

a 8 + 18a 4 + 24a 2 + 21 - 4(1 + a 2 )(l - a 2 fp = (4.1) 



where —3/2 < p < 3 [32|. Because the coordination number of the kagome lattice is even, 
this locus is symmetric under a — > —a. In Ref. fLl|, the locus was plotted in the z and u = z 2 



planes (see also Ref. |4!|). Here it is shown as the gray curves in the a plane, consisting of 
a "dumbell" part and a complex conjugate pair of circular arcs which intersect the dumbell 
at four multiple points (the analytic expressions for which are given in Ref. |ll[|). The inside 
of the dumbell region is the PM phase, the c. c. regions between the narrow neck of the 
dumbell and the circular arcs are O phases, and the region outside of B and extending to 
complex infinity is the (CT extension of the) FM phase. The PM-FM critical point is given 
by a ka g,PM-FM,q=2 = -^kag,t A =2 = 3 1/4 (2 - V3y 1/2 = 2.542... Just there is no physical AFM 
phase (owing to the frustration of the Ising AF on the kagome lattice), so also there is no 
complex-temperature extension thereof. Comparing the zeros calculated with the different 
boundary conditions, we find that with the (pbc,fbc) choice, the zeros on the neck of the 
dumbell and on the outer circular arcs track the exact curves well, while those on the right 
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(left) lie slightly outside (inside) the CT phase boundaries. With the (pbc,fbc) boundary 
conditions, not all sites have even coordination number, so that Z contains some odd powers 
of a, and hence the a — > —a symmetry of the exact boundary is not precisely maintained by 
the zeros. In passing, we note that because the sites on the upper and lower boundaries have 
odd coordination number A = 3, theorem 6 of Ref. Jll| implies that Z(z = —1) = 0, and 
this zero (which is multiple) is evident in Fig. |[ For the (pbc,fbc)A boundary conditions, 
(i) most of the zeros near to the dumbell lie closer to the exact curves, but the zeros near 
the arcs lie farther away from them, as compared with the situation for the (pbc,fbc) choice; 
and (ii) Z, and hence its zeros, is invariant under the negation a — > —a, in contrast to case 
with the (pbc,fbc) case; (hi) because all sites have even coordination number, there is no 
zero in Z at z = — 1. For both types of boundary conditions, the density of zeros in the 
vicinity of the PM-FM critical point ak a g,PM-FM,q=2 decreases in a manner consistent with 
the form from eq. ( |2.6|) with a = for the 2D Ising model, viz., g ~ s as s — > 0, where s the 
arclength along B away from apM-FM,q=2- 



B. q = 3 Case 



Re (a) 

FIG. 11. Zeros of Z in the a plane for the q = 3 Potts model on a kagome lattice of size 4x8 
hexagons and (pbc,fbc) boundary conditions. 
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FIG. 12. Zeros of Z for the q = 3 Potts model on a kagome lattice of size 4x8 hexagons and 
(pbc,fbc)A boundary conditions. 

We show our zeros of Z for the q 



3 Potts model on the kagome lattice in Figs. 11 



and In this case we use lattices of sizes N x x N y with N y larger than N x in order to 
compensate for the fact that the free boundaries are in the y direction and free, as contrasted 
with periodic, boundary conditions introduce greater finite-size effects. The zeros suggest 
that in the thermodynamic limit, the inferred CT phase diagram for the q = 3 kagome lattice 
may involve somewhat simpler boundary curves than was the case for the same model on the 
honeycomb lattice. There is a high-temperature PM phase, a low-temperature FM phase, 
and there are strong indications of a third CT phase whose right-hand boundary crosses the 
real axis at a = 0, corresponding to a zero-temperature critical point of the q = 3 Potts 
antiferromagnet on this lattice. This is in good agreement with the known property that this 
model does have such a T = critical point (which can be related to the T = critical point 
of the q = 4 Potts antiferromagnet on the triangular lattice) [ 4"3"| , |4"4" |. The inferred position 
where the CT boundary crosses the real axis on the right is at ak ag ,PM-FM, q =3 = 2.84(4). This 
is in accord, to within the uncertainty, with the value of ak a g,PM-FM,q=3 = 2.87646(4), i.e., 
Zkag,PM-FM, q =3 = 0.347650(5)) obtained from series analysis ||. The left-hand boundary 
of the third phase crosses the real axis at about a = —0.96(3). From Fig. [T2|, where the 
boundary between the (CT extensions of the) PM and the FM phases is probably best 
represented, we infer that the leftmost point where this CT phase boundary crosses the real 
axis is at ak ag ,e, q =3 = —2.54(6). This point is manifested as a singularity in thermodynamic 
quantities evident in low-temperature series analysis, which obtains ak ag ,e,q=3 = —2.486(3) 
(i.e., Zkag,e,q=3 = —0.4023(5)). Although the CT phase boundary is not symmetric under 
a — > —a as was true for q = 2 on this lattice, one can still discern a remnant of the dumbell 
shape that occurred for the q = 2 case. As before for the honeycomb lattice with q = 3, 
the g ~ s 2 / 3 decrease in the density of zeros in the vicinity of the PM-FM critical point is 



18 



consistent with the calculated zeros. 



C. Case q = 4 



Im(a) 



\ 



Re (a) 

FIG. 13. Zeros of Z for the q = 4 Potts model on a kagome lattice of size 4x8 hexagons and 
(pbc,fbc) boundary conditions. 
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FIG. 14. Zeros of Z for the q = 4 Potts model on a kagome lattice of size 4x8 hexagons and 
(pbc,fbc)A boundary conditions. 



For q = 4, we present our results in the Figs. [13] and |14[ The main differences between 
the locus of zeros, as compared with the case of q = 3 is, first, that the crossing which occurs 
at a = for q = 3 is shifted to a negative value, a = —0.21(2) for q = 4. This demonstrates 
that the q = 4 Potts antiferromagnet on the kagome lattice has no PM-AFM phase transition 
(or any hypothetical non-symmetry breaking phase transition) at finite temperature or any 
critical point at T = 0. This conclusion also follows for q = 5 since increasing q beyond 4 
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has the effect of making the spins "floppier" and the model more disordered. For q > 6, 
this conclusion has been proved rigorously p8| . Second, it appears that the previously 
presumably closed inner ring of zeros has now opened at its leftmost point, which would 
imply that now there would be only two phases (with their CT extensions), the PM and 
FM. The values inferred for ak ag ,pM-FM,q=4 and ak a g,e, q =4 are in accord with the values 
obtained from series analysis Since the specific heat critical exponent a = 2/3 for the 
q = 4 Potts model on 2D lattices, eq. Q2.6|) gives g ~ s 1 ^ 3 for the manner in which the 
density of zeros vanishes as one approaches the PM-FM critical point along the CT phase 
boundary. In particular, this implies that the decrease in density should be less rapid for 
q = 4 than for q = 3, and, indeed, this is evident from a comparison of our plots of zeros for 
these two cases on the kagome lattice. 



D. Partition Function Zeros on the Diced Lattice 

As with the honeycomb lattice and its dual, the triangular lattice, our zeros, in the a 
plane, of the partition function for the g-state Potts model on the kagome lattice also yield 
equivalent zeros of the same model on the lattice that is dual to the kagome lattice, in the 
plane of the variable ad given in eq. (|2.5| ). Henceforth, we shall suppress the subscript d on 
ad- This dual lattice is called the diced lattice; as discussed above, it is not an Archimedean 



lattice ||36|| ; rather, it is a tiling of the plane with identical rhombi such that, as one traverses 
a circuit along the edges of each rhombus, one passes vertices with coordination number 
3,6,3,6 in sequence. Thus, in standard mathematical notation, the diced lattice is the lattice 
[3 ■ 6 ■ 3 ■ 6] dual to the (3 • 6 ■ 3 ■ 6) (= kagome) lattice. Some relevant properties of the 
diced lattice are noted in (Table II of) Ref. |27j. Although the faces of the diced lattice 
are identical, the vertices are not (this is the dual of the property that the vertices of an 
Archimedean lattice are identical but the faces are, in general, not, since an Archimedean 
lattice can consist of more than one type of regular polygon). In particular, the diced lattice 
has vertices of two different types: one with an odd degree (= coordination number) A = 3, 
and the other with even degree, A = 6. Indeed, the diced lattice is bipartite, and its two 
sublattices, which we may denote A 3 and A 6 , are comprised of the vertices with degree A = 3 
and A = 6, respectively. The vertices in the A 3 and A 6 sublattices occupy the respective 
fractions fs = 2/3 and f$ = 1/3 of all the vertices. This is quite different from bipartite 
Archimedean lattices, where the vertices on each of the two sublattices occupy the same 
fraction, / = 1/2, of the total number of vertices (as a consequence of the fact that on an 
Archimedean lattice, all vertices are equivalent). 
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FIG. 15. Zeros of Z for the Ising model on a diced lattice, obtained via duality from a kagome 
lattice of size 4x8 hexagons and (pbc,fbc)A boundary conditions. 

For our discussion of the CT phase diagrams of the g-state Potts model with q = 3,4 
on the diced lattice, it is instructive to begin by discussing the q = 2 case, for which one 
can use exact results on the free energy. The CT phase diagram is shown in Fig. [15j. The 
exact CT phase boundary B [H2J,|TTJ is shown as the dark solid curve. Using duality and the 
z — > —z symmetry of the boundary for the Ising model on the kagome lattice (the latter of 
which follows from the even coordination number of that lattice), it follows that the B in 
Fig. 13 for the Ising model on the diced lattice is the same as B in the v plane for the model 
on the kagome lattice, where v — (1 — z)/ (1 + z) . Note that although the physical FM and 
AFM phases are disjoint, the respective complex-temperature extensions of these phases are 
analytically connected. The reason for this is that, in contrast to bipartite Archimedean 
lattices, the two sublattices A 3 and A 6 of the diced lattice do not occupy the same fraction 
of the total lattice. Thus, reverting to conventional Ising model notation for this discussion, 
let us define M\ 3 and M\ 6 as the magnetizations of the sublattices A 3 and A 6 and M uni f, 
and M st ag. as the uniform and staggered magnetizations, all per unit area of the total lattice, 
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with 



M unif . = M As + M A6 (4.2) 

M stag . = M A3 - M Ae (4.3) 

If the present lattice had been Archimedean, with each sublattice occupying a fraction 1/2 
of the total, then M stag . would vanish identically not just in the PM phase but also in the 
FM phase, and M un if_ would vanish identically not just in the PM phase but also the AFM 
phase, so that the FM and AFM phases, and their complex-temperature extensions, could 
not be analytically connected with each other. However, because the sublattices of the diced 
lattice occupy different fractions of the total lattice, it follows that in the limit of complete 
sublattice spin ordering, M As =2/3 and M Ae = 1/3 and hence, besides the obvious result, 
M uni fXa = oo) = 1, one has 

M um/ .0 = 0) = ^ (4.4) 

M stag .{a = oo) = i (4.5) 

M stag .(a = 0) = 1 (4.6) 

That is, the uniform magnetization M un if. does not vanish even in the region of complete 
sublattice magnetizations of opposite sign, at T = for J < 0, i.e., a = 0, and the staggered 
magnetization M stag . does not vanish even in the limit of complete sublattice magnetizations 
of the same sign, at T = for J > 0, i.e., a = oo. Hence, there exist paths that connect 
the points a = oo and a = in the complex a plane. Of course, if one restricts to the 
physical temperature interval < a < oo, then the physical FM and AFM phases cannot 
be analytically connected, since they are separated by the PM phase, where both M uni ^ 
and M stag , vanish identically. However, the complex-temperature extensions of the FM and 
AFM phases are analytically connected, as is shown by the existence of the paths alluded to 
above. 

The CT phase diagram is thus as follows: (see Fig. ITB: first, there is a symmetric, high- 



temperature PM phase around the point a = 1 that includes the interval a\ < a < a 1 1 on 
the real axis, where 



Cti 

2 



-(1 + V3) \l - (2V3 - 3) 1/2 1 = 0.43542... (4.7) 
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Secondly, there is the single complex-temperature extension of the two different physical FM 
and AFM phases; this extension includes the intervals — oo < a < a\ and a^ 1 < a < oo 
on the real axis (see eq. (4.12) in Ref. [|ll|]) and extends outward to complex infinity in the 
a plane. We label this phase as CT(A)FM. Third, there is a complex-conjugate pair of O 



phases. In Fig. [L5| we have shown the zeros computed with one particular set of boundary 
conditions; in this case and also with the other types of boundary conditions, these zeros 
agree well with the exact results. (This is dual to the same statement for the q = 2 kagome 
lattice.) 

Proceeding to the cases that have not been exactly solved, in Fig. [16] we show our 
zeros for the q = 3 Potts model on the diced lattice, obtained via duality from those on the 
kagome lattice. For this and q = 4, we show results with only one set of boundary conditions, 
since the other boundary conditions yield similar results. As noted in Ref. |fL3"l , from the 
finding in Ref. of a CT singularity at Zt = —0.4023(5) in the q = 3 Potts model on the 
kagome lattice, it follows, using the duality connection, that the q = 3 Potts model has a 
phase transition from the PM phase to the FM-AFM phase at the point adi ce d,q=3,PM-AFM = 
0.1393(8). This constitutes the left border of the physical PM phase on the positive real 
a axis. Moreover, again by duality, from the PM-FM transition point of the model on the 
kagome lattice, determined from series analysis in Ref. || to be at z c = 0.347650(5), it follows 
that the position of the PM-FM transition of the q = 3 Potts model on the diced lattice is 
at adiced,PM-FM, q =3 = 2.59876(4). As in the q = 2 case, and for the same reason, although 
the physical FM and AFM phases are disjunct, their complex-temperature extensions are 
analytically connected. We thus again label this extension as the CT(A)FM phase. The other 
CT phases include the extension of the PM phase and an O phase in the Re(a) < half- 
plane. Our finding that the CT phase boundary for the q = 3 Potts model on the kagome 
lattice has a component that passes through a = 0, corresponding to a zero-temperature 
critical point in that model, implies, by duality, that the boundary of the O phase in the 
model on the diced lattice crosses the real a axis on the left at a = —2. 

We show our zeros for the q = 4 Potts model on the diced lattice in Fig. [17|. For 
this case, from the value ak ag ,PM-FM, q =i = 3.1561(5) obtained from series analysis in 
Ref. 0, we deduce, using duality, that the PM-FM critical point for the diced lattice is 
^diced,PM-FM, q =i = 2.8552(5). Further, from the value obtained for the CT singularity, 
Zkag,q=i,e = —0.42 ± 0.01 J|, we have deduced, again using duality, that the q = 4 Potts 
antiferromagnet on the diced lattice has no finite-T phase transition and also is not critical 
at T = 0, since 



V{a kam= ^) = -(0.18 ±0.02) (4i 
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is negative |13| . There is thus no AFM phase for q = 4. Because increasing the value of q for 
a fixed temperature makes the spins floppier, this result implies that there is also no AFM 
phase for q > 5. In the context of the complex-temperature phase diagram, the point in eq. 
( |4.8|) corresponds to the point where the two arcs close in the thermodynamic limit and the 
left-hand boundary of the PM phase crosses the real a axis in Fig. [17|. In this figure one 
also sees a curve in B in the Re(a) < half-plane. 
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FIG. 16. Zeros of Z for the q = 3 Potts model on a diced lattice with (pbc,fbc) boundary 



conditions and of size equivalent to 4 x 8 hexagons on a kagome lattice. 
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FIG. 17. Zeros of Z for the q = 4 Potts model on a diced lattice with (pbc,fbc) boundary 
conditions and of size equivalent to 4 x 8 hexagons on a kagome lattice. 



V. CONCLUSIONS 

We have calculated complex-temperature zeros of the partition function for the g-state 
Potts model on the honeycomb and kagome lattices. These results give useful information 
about the complex-temperature phase diagrams and singularities of these models. 
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